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Abstract 

Inhomogeneities in the initial QCD matter density distribution increase the production of thermal 
photons significantly compared to a smooth initial-state-averaged profile in the region pj > 1 
GeV/c in an ideal hydrodynamic calculation. This relative enhancement is more pronounced for 
peripheral collisions, for smaller size systems as well as for lower beam energies. A suitably 
normalized ratio of central-to-peripheral yield of thermal photons reduce the uncertainties in the 
hydrodynamical initial conditions and can be a useful parameter to study the density fluctuations 
and their size. The fluctuations in the initial density distribution also lead to a larger elliptic flow 
of thermal photons for pj > 2.0 GeV/c compared to the flow from a smooth profile. 

1. Thermal photons from event-by-event hydrodynamics 

Event-by-event hydrodynamics with fluctuating initial conditions (IC) leads to better agree- 
ment with experimental data for particle spectra and elliptic flow of hadrons produced in heavy 
ion collisions than hydrodynamics with a smooth initial-state-averaged profile QJJ] . Thermal pho- 
tons with pj > 1 GeV/c are mostly emitted from the hot and dense early stage of the system 
expansion. They are thus especially suitable for probing IC fluctuations. For 0-20% central 
Au+Au collisions at RHIC, the thermal photon pj spectrum from fluctuating IC using ideal hy- 
drodynamics is found to be about 10% flatter than the spectrum from smooth IC in the region 
2 < pr < 4 GeV/c. Consequently, it explains the PHENIX direct photon data better in that pj 
range (2J]. We show that a systematic study of photon observables at different collision central- 
hies and collision energies provides useful information about the density fluctuations and their 
size in the IC (3[ . 

The event-by-event hydrodynamic framework developed in IH is used here (as in J2l) to cal- 
culate the thermal photon pj spectra at different collision centralities at RHIC and LHC energies 
as well as photon elliptic flow from fluctuating IC. This is a Monte Carlo Glauber based model 
where the entropy density s is distributed around the wounded nucleons (WN) in the transverse 
plane using a 2D Gaussian function of the form 
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Figure 1: [Left] pj spectra of thermal photons from smooth and fluctuating IC at RHIC for different collision centralities. 
[Right] Relative enhancement due to fluctuations in the IC at RHIC in comparison with LHC for 0-20% centrality bin 
(from (J). 



Here, the free parameter <x determines the fluctuation size. We use state of the art photon rates 
(QGP and hadronic rates from 01 and H respectively) to calculate the photon pj spectra at 
RHIC and LHC energies for different collision centralities. The initial formation time to of the 
plasma is taken as 0.17 and 0.14 fm/c at RHIC and LHC, respectively, from the EKRT model J3|]. 
Results from the fluctuating IC are obtained by taking a final state average of photon pj spectra 
from a sufficiently large number of random events. At different centrality bins the smooth IC 
is constructed by taking an initial-state average of 1000 fluctuating profiles, which essentially 
removes all the fluctuations from the interior of the fireball. 

The pj spectra of thermal photons from Au+Au collisions at RHIC are shown in the left 
panel of FigureQ]for three different centrality bins. The results from the smooth and fluctuating 
IC are obtained using fixed values of cr{- 0.4 fm) and to (= 0.17 fm/c) for all the centrality bins. 
The enhancement in photon production due to IC fluctuations is found to be more pronounced 
for peripheral collisions than for central collisions, i.e., the relative importance of the 'hotspots' 
increases for peripheral collisions. The presence of even a few hotspots in the fluctuating IC 
enhance the photon production significantly as compared to the smooth IC as the photon yield 
falls rapidly towards peripheral collisions. Similar to RHIC, the effect of the IC fluctuations is 
also found to be stronger for peripheral collisions than for central collisions at the LHC. However, 
the relative enhancement is found to be less at the LHC than at RHIC for the same centrality bin. 
The ratio of the thermal photon yield from the fluctuating and smooth IC at RHIC and LHC 
energies for the 0-20% centrality bin is shown in the right panel of Figure Q] 

1.1. Initial formation time dependence 

The thermal emission of photons depends strongly on the initial formation time of the plasma. 
A smaller to leads to a larger initial temperature and more high pj photons from the plasma 
phase. However, the value of To is not known precisely and it varies within a range of 0.17 to 
0.60 fm/c in different hydrodynamic calculations at RHIC energy. We see that the pj spectra 
from smooth as well as from fluctuating IC fall sharply in the range 2 < pj < 4 GeV/c for all 
centrality bins when to is changed from 0.17 fm to 0.60 fm/c (keeping the total entropy of the 
system fixed, see Figure 4 in [3]). In addition, the value of to may not remain the same for 
different centrality bins. The system is more dilute for peripheral collisions and one can expect a 
larger To there compared to central collisions. 
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Figure 2: [Left] pj spectra are sensitive to the value of To which may also vary with collision centrality. [Right] Ratio of 
central to peripheral yield of thermal photons as a function of collision centrality (from (3j). 



The left panel of Figure[2]shows pj spectra of thermal photons at the LHC energy considering 
centrality dependent To values (calculated from a non-local version of the EKRT model JfJ with 
^nlo° = anc ' fi xe d T o f° r ai l centrality bins (see yfl for details). A larger to for peripheral 
collisions decreases the production compared to a smaller to, whereas IC fluctuations increase 
the production. However, due to the uncertainties in the initial formation time it is difficult to 
conclude anything about the density fluctuations and their size in the IC by looking at the pj 
spectra alone. 

1.2. Ratio of central to peripheral yield R y cp 

A suitably normalized ratio of central to peripheral yield can be a useful measure of the 
fluctuation size parameter by reducing the uncertainties in the initial conditions. We define a 
quantity Rl p as yfl 

r7 | _ dN/d 2 p T dY\ Q - m% ^ N bi „\i-j% ^ 
cpl dN/d 2 p T dY\ H% N bin \ -io%' 

where the value of i is changed from 10 to 70 in steps of 10 and j-i+10 in the following. For 
pQCD direct photons (not included here) such a ratio is close to unity. A^in is the number of 
binary collisions for a particular centrality bin. 

The right panel of Figure [2] shows the variation of Rl p as a function of collision centrality 
for 200A GeV Au+Au collisions at RHIC and for different values of pj and <x. For smaller pj 
and for central collisions the value of R cp is close to unity by definition. For pj > 2 GeV/c, RL 
rises rapidly with collision centralities and a significant difference between the results from the 
smooth and fluctuating IC can be observed even for larger values of <x and To JH. 

1.3. Elliptic flow from fluctuating initial conditions 

The elliptic flow of direct photons is believed to be dominated by the flow contribution from 
thermal radiation up to pj ~ 5 GeV/c as the contributions from all other sources (apart from 
thermal) are expected to be small in that pj range. We calculate the elliptic flow of thermal pho- 
tons from event-by-event hydrodynamics with respect to the participant plane using the relation 
v>2 = (cos 2(<p - 0pp))event> where (ppp is the participant plane angle [7]. The obtained \'2{pt) from 
smooth and fluctuating IC at RHIC for the 20^-0% centrality bin is shown in the left panel of 
Figure [3] Fluctuations in the initial density distribution lead to a larger transverse flow velocity 

3 



0.(16 
0.05 
_ 0.04 
0.03 
0.02 
0.111 



Elliptic Bow of thermal photons 

U00A GeV Au+Au(aRHIC 
20-40% Centrality bin 
<j=0.4 fm, t„= 0.17 fm/c 



~ v 2 (PP) 

— v 2 (RP) 

— smooth (MC) 



i Thermai Photons " ' 

z 200A GeV Au+Au@RHIC 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 

p T (GeV/c) 



p T (GeV/c) 



Figure 3: [Left] Elliptic flow of thermal photons at RHIC from fluctuating initial conditions. [Right] Time evolution of 
thermal photon pr spectra and elliptic flow (normalized by the final pj spectra and elliptic flow respectively). 



and larger elliptic flow especially for pj > 2 GeV/c compared with a smooth profile. However 
the results from the fluctuating IC are still found to be well below the PHENIX data QSJ] - The 
right panel of Figure [3] shows the time evolution of the thermal photon pj spectra and elliptic 
flow (normalized by the final pj spectra and elliptic flow respectively) of a single event. We see 
that at pj — 4 GeV/c, more that 80% of the photons are emitted before first one fm, however less 
than 10% of elliptic flow is developed by that time. 

2. Summary and conclusions 

In conclusion, the thermal photon pj spectra from smooth and fluctuating IC are compared 
for different collision centralities at the RHIC and LHC energies. The relative enhancement due 
to fluctuations in the IC is found to be more pronounced for peripheral collisions and for lower 
beam energies. A suitably normalized ratio of central-to-peripheral yield as a function of pj 
and collision centrality can be a useful measure of the fluctuation size scale by reducing the 
uncertainties in the initial condition of a hydrodynamic calculation. Fluctuations in the IC also 
increase the elliptic flow significantly compared to a smooth profile for pj > 2 GeV/c. 
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